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Abstract 

Background: Mycobacterium abscessus group includes antibiotic-resistant, opportunistic mycobacteria that are 
responsible for sporadic cases and outbreaks of cutaneous, pulmonary and disseminated infections. However, 
because of their close genetic relationships, accurate discrimination between the various strains of these 
mycobacteria remains difficult. In this report, we describe the development of a multispacer sequence typing (MST) 
analysis for the simultaneous identification and typing of M. abscessus mycobacteria. We also compared MST with 
the reference multilocus sequence analysis (MLSA) typing method. 

Results: Based on the M. abscessus CIP104536 1 genome, eight intergenic spacers were selected, PCR amplified and 
sequenced in 21 M. abscessus isolates and analysed in 48 available M. abscessus genomes. MST and MLSA grouped 
37 M. abscessus organisms into 12 and nine types, respectively; four formerly "M. bolletii" organisms and M. abscessus 
Ml 39 into three and four types, respectively; and 27 formerly "M. massiliense" organisms grouped into nine and five 
types, respectively. The Hunter-Gaston index was off 0.912 for MST and of 0.903 for MLSA. The MST-derived tree 
was similar to that based on MLSA and rpoB gene sequencing and yielded three main clusters comprising each the 
type strain of the respective M. abscessus sub-species. Two isolates exhibited discordant MLSA- and rpoB gene 
sequence-derived position, one isolate exhibited discordant MST- and rpoB gene sequence-derived position and 
one isolate exhibited discordant MST- and MLSA-derived position. MST spacer n°2 sequencing alone allowed for the 
accurate identification of the different isolates at the sub-species level. 

Conclusions: MST is a new sequencing-based approach for both identifying and genotyping M. abscessus mycobacteria 
that clearly differentiates formerly "M. massiliense" organisms from other M, abscessus subsp. bolletii organisms. 

Keywords: Mycobacterium, Mycobacterium abscessus, Mycobacterium massiliense, Mycobacterium bolletii, Multispacer 
sequence typing, Genotyping 



Background 

Mycobacterium abscessus mycobacteria are increasingly 
being cultured from respiratory tract specimens col- 
lected from patients with chronic pulmonary diseases, 
including cystic fibrosis [1-9]. These mycobacteria are 
also responsible for skin and soft-tissue infections 
following surgical and cosmetic practices [10-12] and 
catheter-related bacteremia [13,14]. These infections are 
particularly critical for immune-compromised patients 
and may be fatal [15]. Water is suspected as a source 
of infection, as M. abscessus mycobacteria have been 
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isolated from tap water [16]. Moreover, M. abscessus 
mycobacteria have been shown to be resistant to water- 
borne free-living amoebae [17,18]. M. abscessus infec- 
tions are also associated with treatment failure owing, 
due to the natural broad-spectrum resistance to antibio- 
tics in addition to acquired resistance, with subtle differ- 
ences in the antibiotic susceptibility pattern being 
observed among isolates [19]. 

Indeed, M. abscessus is comprised of a heterogeneous 
group of mycobacteria currently classified into M. 
abscessus subsp. abscessus and M. abscessus subsp. bolle- 
tii [20,21], with the later subspecies accommodating 
mycobacteria previously identified as "Mycobacterium 
bolletii" or "Mycobacterium massiliense" [18,22]. How- 
ever, these organisms are nearly indistinguishable using 
phenotypic tests including the mycolic acid pattern 
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analysis and share 100% 16S rRNA gene sequence similar- 
ity [20]. They were initially differentiated on the basis of 
>3% rpoB gene sequence divergence and different anti- 
microbial susceptibility patterns [23,24]. Nevertheless, 
confusing results based on rpoB sequencing have been 
reported [21], and combining sequencing of the rpoB, 
hsp65 and secA genes has been advocated for the optimal 
identification of the M. abscessus mycobacteria [25] . 

To further decrypt the diversity and genetic relation- 
ships among M. abscessus organisms, we investigated a 
collection of reference, sequenced genomes and clinical 
M. abscessus isolates using multispacer sequence typing 
(MST), which is a sequencing-based approach previously 
used for the species identification and genotyping of 
Mycobacteria, including Mycobacterium avium [26] and 
Mycobacterium tuberculosis [27] and non-mycobacterial 
pathogens, such as Yersinia pestis [28], Rickettsia prowa- 
zekii [29] and Bartonella quintana [30]. This approach 
was here compared with multilocus sequence analaysis 
which relies the sequencing of 5-8 genes (21, 25), and 
rpoB genes sequencing (23, 24). 

Methods 

Bacterial isolates 

Reference M. abscessus CIP104536 T , M. abscessus 
DSMZ44567 (German Collection of Microorganisms 
and Cell Cultures, Braunschweig, Germany), M. absces- 
sus subsp. bolletii CIP108541 T (herein referred as 
"M. bolletii") and M. abscessus subsp. bolletii CIP108297 T 
(herein referred as "M. massiliense" [23]) were used in this 
study. In addition, a collection of 17M. abscessus clinical 
isolates from the mycobacteria reference laboratory of the 
Mediterranee Infection Institute, Marseille, France were 
also studied (Table 1). All of the mycobacteria were grown 
in 7H9 broth (Difco, Bordeaux, France) enriched with 10% 
OADC (oleic acid, bovine serum albumin, dextrose and 
catalase) at 37°C. As for the identification, DNA extraction 
and rpoB partial sequence-based identification were per- 
formed using the primers MYCOF and MYCOR2 (Table 1) 
as previously described [24]. In addition, the rpoB gene 
sequence retrieved from 48 M. abscessus sequenced 
genomes was also analysed (Additional file 1) (http:// 
www.ncbi.nlm.nih.gov/) . 

Reference MLSA typing 

Fragments from five housekeeping genes argH (arginino- 
succinate lyase), cya (adenylate cyclase), murC (UDP 
N-acetylmuramate-L-Ala ligase, pta (phosphate acetyl- 
transferase) and purH (phoshoribosylminoimiazolcarboxy- 
lase ATPase subunit) were amplified using the sets of 
primers as previously described (21). The sequences of 
each one of these five housekeeping genes retrieved from 
48 M. abscessus sequenced genomes, were also included 
in the MLSA analysis (Additional file 1). 



MST analysis 

Sequences of the whole intergenic spacers were 
extracted from the reference M. abscessus CIP104536 T 
(ATCC19977) genome (GenBank accession CU458896.1) 
using the perl script software and a total of 8 spacers 
with a 200-700-bp sequence size were further used in 
analysis. For each of these 8 spacers, specific PCR pri- 
mers were designed using Primer3 software v 0.4.0 
(http://frodo.wi.mit.edu/primer3) and tested in silico for 
specificity using BLAST software (http://www.ncbi.nlm. 
nih.gov). The PCR conditions were first optimized using 
DNA extracted from the reference M. abscessus, "M. 
bolletii" and "M. massiliense" isolates before analysis of 
DNA extracted from the 17 clinical isolates (Table 1). 
The PCR amplifications were performed in a 50 ul PCR 
mixture containing 5 ul lOx buffer (Qiagen, Courtaboeuf, 
France), 200 mM each dNTP, 1.5 mM MgCl 2 , 1.25 U Hot- 
Star Taq polymerase (Qiagen), 1 ul each primer (10 pM), 
33 ul nuclease-free water and 5 ul DNA template. The 
amplification program consisted of an initial 15 min 
denaturation step at 95°C followed by 40 cycles of 30 s at 
95°C, 30 s at 60°C and 1 min at 72°C; the amplification 
was completed by a final 5-min elongation step at 72°C 
Negative controls consisting of PCR mixture without 
DNA template were included in each PCR run. The pro- 
ducts were visualized by gel electrophoresis, purified using 
a MultiScreen PCR filter plate (Millipore, Molsheim, 
France) and sequenced in both directions using the 
BigDye Terminator sequencing kit (Applied Biosystems, 
Villebon-sur-Yvette, France), as previously described [27]. 
The sequences were edited using the ChromasPro soft- 
ware (version 1.42; Technelysium Pty Ltd), aligned using 
Clustal W (MEGA 5 software) and compared with the 
reference M. abscessus ATCC 19977 sequences (GenBank 
accession CU458896.1). For MST and MLSA discrimin- 
ation power was calculated using the Hunter-Gaston 
Index [31]: 



DI= 1 



1 



N(N ■ 



where D is the numerical index of discrimination, N is the 
total number of isolates in the sample population, s is the 
total number of different types and nj is the number of 
isolates belonging to the ;'th type. 

Phylogenetic analysis 

Phylogenetic trees were constructed based on rpoB gene, 
concatenated MLSA genes, concatenated spacers and 
MST spacer n°2 sequences using the neighbor-joining 
method with Kimura's two-parameter (K2P) distance 
correction model with 1000 bootstrap replications in the 
MEGA version 5 software package [32]. The rpoB gene 
sequence-based tree was rooted using M. chelonae strain 
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Table 1 Spacers characteristics used in this study 



Name 


Genome position* 


Framing genes* 


PCR primers 


PCR product size (bp) 


Spacer 1 


106145-106396 


MAB_01 04:enoyl-CoA hydratase/isomerise 


F : GGGATGCGCAGATGACGGGG 


506 






MAB_01 05c:oxidoreductase 


R : GCTACCCCGAATGGGGCACG 




Spacer 2 


173727-173985 


MAB_01 76:antigen 85-A precursor 


F : TCGAGTTTCCTCCGGGCGGT 


438 






MAB_0177:antigen 85-A/B/C precursor 


R: AATCCAGGCAGAACGGCCGC 




Spacer 3 


422777-423027 


MAB_0423c:hypothetical protein 


F: GCCATTGCTGTCCGTGCGGT 


344 






MAB_0424:putative protease 


R : GCCGCGAACAGGCCAAACAG 




Spacer 4 


49441 1-494670 


MAB_0495c:hypothetical protein 


F: CGCCCTTGCGCAGGAGTGAT 


528 






MAB_0496c:hypothetical protein 


R: GCCTGGTTCGGACGGTGACG 




Spacer 5 


761805-762060 


MAB_0761c:putative 3-hydroxyacyl-CoA dehydrogenase 


F : ACCACATCGGCGAGCGTGTG 


545 






MAB_0762:hypothetical protein 


R : CCAACACCGGGTCGCGGTAC 




Spacer 6 


771 170-771436 


MAB_0772c:hypothetical protein 


F : CGTCGGTCTTGCCGACCGTC 


600 






MAB_0773:hypothetical protein 


R : GGCGCCGACGATCTAGCACC 




Spacer 7 


880381-880639 


MAB_0887c:hypothetical protein 


F: CGGCAGTGCAAGGTGCGTTG 


519 






MAB_0888c:putative fumarylacetoacetase 


R : GCACCGTGTCCGGTCCTCAG 




Spacer 8 


959422-959678 


MAB_0950c:putative amino acid permease family protein 


F: GGGGCGTATGCGCCGTTACC 


474 






MAB_0951:putative aminoglycoside phosphotransferase 


R : CGAACGCGCTGTGATTCGGC 




Spacer 9 


1002935-1003200 


MAB_0995:hypothetical protein 


F : GGCCGCGACAAGCTGATCGT 


684 






MAB_0997c:hypothetical protein 


R: ATGCAGGGCACCGTGCGTAG 




Spacer 10 


1216613-1216879 


MAB_1201c:transcription elongation factor GreA 


F: CGTTCTCGCGCAGGTCTCCC 


517 






MAB_1202c:hypothetical protein 


R: CCGAACGATCCGTGCCGGTC 




Spacer 1 1 


1818877-1819188 


MAB_1818:hypothetical protein 


F: AGCCAACTGCCATGGCGCTT 


495 






MAB_1819c:hypothetical protein 


R : ACCGAGACGTCATGCACCGC 





* With reference to M. abscessus ATCC 19977 genome. 

CIP 104535 T and M. immunogenum strain CIP 106684 T 
rpoB gene sequences. A heatmap was constructed using 
the R statistical software based on the spacer profile as a 
distance matrix. 



Results and discussion 

rpoB identification and rpoB tree 

The identification of M. abscessus CIP104536 T , M. 
abscessus DSMZ44567, M. bolletii CIP108541 T and M. 
massiliense CIP108297 T was confirmed by partial rpoB 
sequencing. The sequences were deposited in the 
GenBank database (GenBank accession: KC352778 - 
KC352795). Isolates PI, P2.1, P2.2, P2.3, P2.4, P2.5, P3.1, 
P3.2, P4, P5, P6, P7 and P8 exhibited 99% rpoB sequence 
similarity with M. abscessus ATCC19977 T and were 
identified as M. abscessus. Isolates P9 and P10 exhibited 
99% rpoB sequence similarity with "M. bolletii" 
CIP108541 1 and were identified as "M. bolletii" whereas 
isolate Pll exhibited 99% rpoB sequence similarity with 
"M. massiliense" CIP108297 T and was identified as "M. 
massiliense". A total of 23 M. abscessus sequenced gen- 
omes were identified as M. abscessus since they exhib- 
ited 98 to 100% similarity with the M. abscessus type 
strain rpoB partial gene sequence. M. abscessus M24 
shared 99% similarity with the M. bolletii type strain 



partial rpoB gene sequence. A total of 26 M. abscessus 
and "M. massiliense" sequenced genomes shared 99% to 
100% similarity with "M. massiliense" partial rpoB gene 
sequence. The tree built from 69 partial rpoB gene 
sequences showed three distinct groups, each compris- 
ing the type strain (Figure la). 



Reference MLSA analysis 

Fragments for the expected size were amplified and 
sequenced for the five MLSA genes. The sequences were 
deposited in the GenBank database (GenBank accession: 
KC352742 - KC352759, KC352760 - KC352777, 
KC352796 - KC352813, KC352814 - KC352831, 
KC352832 - KC352849). Concatenation of the five 
sequences yielded a total of 19 different types, including 
9 types for 37 M. abscessus organisms, four types for 4 
"M. bolletii" organisms and M. abscessus M139 and five 
types for 27 "M. massiliense" organisms. The Hunter- 
Gaston Index for MLSA was of 0.903. The MLSA tree 
based on the five gene concatened sequences showed 
three principal clusters, i.e. a M. abscessus cluster, 
a "M. bolletii" cluster and a "M. massiliense" cluster 
(Figure lb). Latter cluster comprised of five sub-clusters 
with "M. massiliense" type strain and Pll strain sub- 
clustering together close to M. abscessus 5S strain. Also, 
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Figure 1 Phylogenetic tree based on rpoB gene sequence (a); based on the concatenated five MLSA gene sequences (b); and based on 
the concatenated eight polymorphic spacers (c). 



MLSA-derived tree clustered M. abscessus M139 strain 
and P5 strain respectively identified as "M. massiliense", 
close to the "M. bolletii" whereas both strains clus- 
tered with M. abscessus in the rpoB gene sequence- 
derived tree. 

MST analysis 

Analysis of the reference M. abscessus ATCC 19977 
complete genome sequence yielded 3538 intergenic 
spacers with > 300 spacers were 200-700 bp in length. 
Successful PCR sequencing was achieved for 8 spacers in 
all the isolates studied; the sequences were deposited in 
the GenBank database (GenBank accession: KC352850 - 
KC352890). In M. abscessus isolates, including the 37 
sequenced genomes, the spacer sequence variability was 
generated by one to 12 single nucleotide polymorphisms 
(SNPs) (spacers n°l and n"8), one to 18 SNPs and one to 
two nucleotide deletions (spacer n°2), one to two SNPs 
(spacers n°3 and n°7) and nucleotide insertion (spacers 



n°2 and n°5). In "M. bolletii" isolates, the spacer sequence 
polymorphisms were generated by one SNP for spacer 
n°l, two SNPs and one deletion for spacer n°2, two SNPs 
for spacer n°3 and nine SNPs for spacer n°7. In "M. mas- 
siliense" isolates, including 28 sequenced genomes, the 
spacer sequence polymorphism were generated by nine 
SNPs and one insertion (spacer n°l), one insertion (spacer 
n°3), five SNPs and two insertions (spacer n°4), one SNP 
(spacer n°5) and two SNPs (spacer n°7). Concatenation of 
the eight spacer sequences yielded a total of 24 types, 
with the 37 M. abscessus organisms grouped into 12 spa- 
cer types, four formerly "M. bolletii" organisms grouped 
into three spacer types and 28 formerly "M. massiliense" 
organisms grouped into nine spacer types. This yielded a 
Hunger-Gaston Index of 0.912. Spacer n°5 was found to 
be the most variable of the eight spacers under study, 
exhibiting 13 different alleles (Table 2). When combining 
the eight spacer sequences, a unique MST profile for each 
reference isolate was obtained, i.e., MST1 and MST2 for 
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Table 2 Spacers allelic polymorphism and MST a genotypes of M. abscessus, 


"M. bolletii" and 


"M. massiliense" 


isolates 


Isolates 


Spacerl Spacer2 Spacer3 


Spacer4 Spacer5 Spacer6 


Spacer7 


Spacer8 


Genotype 


/W.abscessus_ATCC1 9977_CIP1 04536T 


1 1 1 


1 


1 1 


1 


1 


1 


Mabscessus_DSMZ44567 


2 1 2 


2 


2 1 


2 


1 


2 


P1 


2 1 2 


2 


2 1 


2 


1 


2 


P2.1 


1 2 1 


3 


1 1 


2 


2 


3 


P2.2 


1 2 1 


3 


1 1 


2 


2 


3 


P2.3 


1 1 1 


1 


1 1 


1 


1 


1 


P2.4 


1 1 1 


1 


1 1 


1 


1 


1 


P2.5 


1 1 1 




1 1 


1 




1 


P2.6 


1 1 1 




1 1 


1 




1 


P3.1 


3 1 2 




1 1 


2 




4 


P3.2 


3 1 2 




1 1 


2 




4 


P4 


1 1 1 




1 1 


1 




5 


P5 


1 1 1 




3 1 


2 




6 


P6 


1 1 1 




1 1 


1 




1 


P7 


4 1 2 


4 


4 1 


2 




7 


P8 


4 1 2 


4 


4 1 


3 




8 


M.abscessus_3A-0930-RJA_0930_R 11111111 1 


M.abscessus_3A-0930-S_3A_0930_S 11111111 1 


M.abscessus_3A-0122-S_3A_0122_S 11111111 1 


M.abscessus_3A-0731_3A^0731 11111111 1 


M.abscessus_3A-0122-R_3A_0122_R 11111111 1 


M.abscessus_3A-0119-R_3A_0119_R 11111111 1 


M.abscessus_6G-0728-R_M6G_0728_R 11111111 1 


M.abscessus_6G-0212_M6G_0212 11111111 1 


M.abscessus_6G-l108_6G_1108 11111111 1 


M.abscessus _6G-0728-S_6G_0728_S 11111111 1 


M.abscessus_6G-0125-R_6G_0l25_R 11111111 1 


M.abscessus_6G-0125-S_6G_0125_S 11111111 1 


M.abscessus_4S-01 16-S_4S_01 16_S 


5 1 2 


5 


5 2 


2 


2 


9 


M.abscessus^4S-01 16-R_4S_01 16_R 


5 1 2 


5 


5 2 


2 


2 


9 


M.abscessus_4S-0206_M4S_0206 


5 1 2 


5 


5 2 


2 


2 


9 


M.abscessus„4S-0726-RB^4S„0726_RB 


5 1 2 


5 


5 2 


2 


2 


9 


M.abscessus_4S-0303_4S_0303 


5 1 2 


5 


5 2 


2 


2 


9 


M.abscessus_4S-0726-RA_4S_0726_RA 


5 1 2 


5 


5 2 


2 


2 


9 


M.abscessus_M93 


3 1 2 


6 


6 1 


2 


3 


10 


M.abscessus_M94 


2 1 2 




7 1 


4 


2 


11 


M.abscessus_M1 52 


2 1 2 


7 


7 1 


2 


3 


12 


M6o//eff/_CIP1 08541 T 


6 3 3 


3 


8 1 


5 


2 


13 


P9 


6 3 3 


3 


8 1 


5 


2 


13 


P10 


7 4 1 


3 


8 1 


2 


2 


14 


M.abscessus_M24 


8 3 4 


8 


8 1 


2 


2 


15 


M.massillienseJLW] 08297T 


5 5 5 


9 


9 1 


6 


3 


16 


P11 


5 5 5 


9 


9 1 


6 


3 


16 


M.massiliense _2B-091 2-S_2B_09 12_S 


9 5 6 


10 


10 2 


7 


3 


17 


M.massiliense_2B-030_ M2B_0307 


9 5 6 


10 


10 2 


7 


3 


17 


M.massiliense_2B-0912-R_2B_0912_R 


9 5 6 


10 


10 2 


7 


3 


17 


M.massiliense_2B-0626_M2B_0626 


9 5 6 


10 


10 2 


7 


3 


17 
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Table 2 Spacers allelic polymorphism and MST a genotypes of M. abscessus, "M. bolletii" and "M. massiliense" isolates 

(Continued) 



M.massiliense_2B-123 1_M2B_ 123 


; 




9 


5 


6 


10 


10 


2 


7 


3 


17 


M.massiliense_2B-0 1 07_M2B_0 1 07 




9 


5 


5 


1 0 


1 0 


2 


7 


3 


17 


IVI.lliUbbllltzltbc __ 1 j- 1 D*t-Uj I U_IVI 1 j_ 


1 

/ j'f 


0310 


9 


5 


5 


1 0 


1 0 


2 


7 


3 


17 


M mnccilionco 7C 7^9 HQ1A /1/I7C 


7 C? 


0914 


9 


5 


6 


10 


10 


2 


7 


3 


17 


ivi.massiiiense_ i 3- i dj-uv / j_m / 3_ 




0915 


9 


5 


6 


10 


10 


2 


7 


3 


17 


ivi.massiiiense_ 13-/31 -uvju_ivi 1 j_ 




0930 


9 


5 


6 


10 


10 


2 


7 


3 


17 


HA mstcctlianco HA 1 Q 

ivi.rnubbiiicnsc _ivi i o 






9 


5 


5 


1 0 


1 0 


2 


7 


3 


17 


M.dbscGSSus_Ml 59 






9 


ft 






1 0 


3 


7 


4 


18 


IVI.UUbLcbbUb_ 4 t/ JZO 






9 


5 


6 


6 


1 1 


4 


7 


3 


19 


M.obscessus_M1 72 






1 0 


7 


2 


9 


1 2 


3 


8 


5 


20 


M.abscessu$_M1 54 






10 


7 


2 


g 


12 


3 


8 


5 


20 


M.abscessus_5S- 1215_5S_ 1215 








5 


2 




1 3 


2 


5 


2 


21 


M.abscessus_5S- 1212_5S_ 1212 








5 


2 




1 3 


2 


5 


2 


21 


M.abscessus_5S-08 1 7_5S_08 1 7 








5 


2 


6 


13 


2 


6 


2 


21 


M.abscessus_5S-0708_5S_0708 








5 


2 


6 


13 


2 


6 


2 


21 


M.abscessus_5S-0422_5S_0422 








5 


2 


6 


13 


2 


6 


2 


21 


M.abscessus_5S-0304_5S_0304 








5 


2 


6 


13 


2 


6 


2 


21 


M.abscessusJS-042 1_5S_042 1 








5 


2 


6 


13 


2 


6 


2 


21 


M.abscessus_M1 56 






10 


7 


2 


11 


12 


3 


9 


5 


22 


M.abscessus_M148 






10 


7 


2 


11 


12 


3 


9 


5 


23 


M.abscessus_M139 






10 


5 


2 


11 


14 


3 


10 


3 


24 
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M. abscessus CIP104536 T and M. abscessus DSMZ44567 
respectively, MST13 for "M. bolletii" CIP108541 T and 
MST16 for "M. massiliense" CIP108297 T . At the sequence 
level, we found that MST1 and MST2 genotypes differ by 
at most nine SNPs, whereas MST1 differed from MST13 
by up to 18 SNPs, one insertion and two deletions and 
from MST16 by 14 SNPs, 11 deletions and two insertions 
(supplementary material). The 17 clinical M. abscessus 
isolates were grouped into eight MST types, named 
MST1 to MST8, with five M. abscessus isolates exhib- 
iting the M. abscessus CIP104536 T MST1 genotype and 
one isolate (PI strain) exhibiting the M. abscessus 
DSMZ44567 MST2 genotype. The P9 "M. bolletii" clinical 
isolate yielded the MST 13 genotype in common with the 
reference "M. bolletii" CIP108541 T , whereas the P10 "M. 
bolletii" clinical isolate yielded a unique MST14 genotype 
that differ from MST13 by two SNPs in spacer n°l. M. 
abscessus M24 yielded the MST15 and differed from 
MST13 by four polymorphic spacers. In "M. massiliense" 
nine different profiles were generated MST 16 to MST24. 
The Pll "M. massiliense" clinical isolate shared the 
MST 16 genotype with the reference "M. massiliense" 
CIP108297 . "M. massiliense" 2B isolate, "M. massiliense" 
IS isolate and "M. massiliense" M18 isolate shared the 



same MST profile (MST17). M, abscessus 5S isolate exhib- 
ited the MST21 profile. 

MST based tree and comparaison with rpoB identification 
and MLSA analysis 

The MST-phylogenetic tree clustered isolates from 
patients PI to P8 with M. abscessus reference strain, iso- 
lates from P9 and P10 with "M. bolletii" and isolate from 
Pll with "M. massiliense", in agreement with their rpoB 
sequence-based identification and MLSA analysis 
(Figure lc). The MST, MLSA and rpoB phylogenetic 
trees separated the M. abscessus isolates into three prin- 
cipal clusters depicted by M. abscessus, "M. bolletii" and 
"M. massiliense" isolates (Figure la, b and c). However, 
MST resolved "M. bolletii" cluster into two sub-clusters 
formed by isolate P5 and all of the other M. bolletii iso- 
lates with a 76% bootstrap value, wich is discordant with 
MLSA and rpoB based tree. Each cluster or sub-cluster 
of the M. abscessus isolates corresponded to different 
genotypes. The "M. massiliense" cluster was more dis- 
perse and divided into six sub-clusters with isolate Pll 
and "M. massiliense" type strain sub-clustering alone. 
The results of this analysis were consistent for 67 iso- 
lates and inconsistent for two isolates P5 and M. 
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Figure 2 Heatmap and clustering of M. abscessus mycobacteria under study based in difference of profile. 



abscessus M139. A heatmap incorporating all spacer pat- 
terns into a matrix further demonstrated that spacer n°2 
was the most discriminating spacer (Figure 2). Hence, 
the tree based on the spacer n°2 sequence also discrimi- 
nated the three M. abscessus, "M. bolletii" and "M. mas- 
siliense" clusters (Figure 3). This discrimination potential 
makes spacer n°2 a useful new tool for the accurate 
identification of M. abscessus subspecies. Furthermore, 
these data indicated that it was readily possible to dis- 
criminate isolates that would have been identified as 
"M. bolletii" [26] or "M. massiliense" [23] using a 



previous taxonomy proposal and are now grouped as 
M. abscessus subsp. bolletii according to a recent 
taxonomy proposal [20,21]. 

Conclusion 

We herein developed a sequencing-based MST genotyping 
technique that allows the accurate identification and 
discrimination of M. abscessus mycobacteria. Therefore, 
MST could be added to the panel of molecular methods 
currently available for genotyping M. abscessus mycobac- 
teria, with the advantages that MST is a PCR and 
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sequencing-based technique, thereby providing a robust 
and accurate result without requiring a high DNA concen- 
tration and purity, as is the case for pulsed-field gel 
electrophoresis (PFGE) [5] and randomly amplified poly- 
morphic DNA (RAPD) [33]. Furthermore, MST targets 
intergenic spacers, which undergo less evolutionary pres- 
sure and are thus more variable than the housekeeping 
genes targeted in multilocus sequence typing [21]. Also, 
MST incorporating sequencing is an open approach to 
described new genotypes more versatile than counting the 
number of tandem repeats [34]. We propose that MST 
could be incorporated into a polyphasic molecular 
approach to resolve the phylogenetic relationships of 
difficult-to-identify M. abscessus isolates [35] . Combining 
MST data with phylogenetic analyses clearly indicated that 
M. abscessus heterogeneity spans beyond the current two 
M. abscessus subspecies, as two "M. massiliense" isolates 
were readily discriminated from the other "M. bolletii" 
isolates [21]. These data, therefore, question the current 
nomenclature of M. abscessus mycobacteria, which 
incorporates mycobacteria previously recognized as "M. 
bolletii" and "M. massiliense" as "M. abscessus subsp. 
bolletii". The data presented here indicate that this no- 
menclature masks the underlying diversity of M. 
abscessus mycobacteria, potentially hampering the 
recognition of microbiological, epidemiological and 
clinical particularities that are linked to each subspe- 
cies. The elevation of "M. massiliense" as a new M. 
abscessus subspecies would accommodate the data 
produced in the present study [24]. 

Additional file 



Additional file 1: rpoB and MLSA genes accession Number of 49 
sequenced genomes. 
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